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We study the velocity distribution function for inelastic Maxwell models, characterized by a 
Boltzmann equation with constant collision rate, independent of the energy of the colliding particles. 
• By means of a nonlinear analysis of the Boltzmann equation, we find that the velocity distribution 

^ ' function decays algebraically for large velocities, with exponents that are analytically calculated. 



B 

I 



o 



I. INTRODUCTION 



, Velocity distributions have over-populated high energy tails in many particle systems with in-elastic interactions, 
as has been discovered theoretically, and later observed in laboratory experiments with granular materials |^ . Instead 
of Gaussian tail distributions, many kinetic equations typically predict stretched exponentials, like exp[— in 
driven inelastic hard sphere systems (IHS) ||^, or even tails with a higher overpopulation, like exp[— in the 
' freely evolving IHS fluids 1^,^. These theoretical predictions have been extensively verified in Direct Monte Carlo 
^ . Simulations (DSMC) of the Boltzmann equation |^,^, and the stretched exponentials have been observed in laboratory 
O ' experiments with granular matter on vibrating plates |jl|]. However, the presence and also the absence of over-populated 
, O , , tails depends very strongly on the energy input or on how the system is thermo-statted [^,^ . 

Recently classes of simplified kinetic models have been studied, so-called Maxwell models, which are characterized by 
T-H , a Boltzmann equation with a collision rate that is independent of the relative kinetic energy of the colliding particles. 
^ ' Maxwell molecules are for kinetic theory, what harmonic oscillators are for quantum mechanics, and dumb-bells for 
[ polymer physics. 

Ben-Naim and Krapivsky Q have presented arguments about the non-existence of scaling solutions of the Boltzmann 
equation for the one-dimensional Maxwell model, and argue in favor of multi-time scales However, Puglisi et al. 
T— I have found an exact scaHng solution f{v,t) ~ {l/vQ{t))f{v/vo{t)) for that equation, i.e. /(c) = (2/7r)(l -|-c^)^. 

T— I ' It does have a power law tail (1/c^). Puglisi et al. also solve the spatially homogeneous Maxwell-Boltzmann 
. equation using MC simulations in one- and two- dimensions, and, more importantly, they show that an arbitrary 
' initial distribution approaches this scaling solution, with power law tails at high energies. In two-dimensions the 
^ , exponent a of the power law tails depends on the degree of in-elasticity, i.e. on the coefficient of restitution a. The 
^ ' goal of this article is to derive these power laws from the dominant small-A: singularity in the Fourier transform of the 
^ ^ velocity distribution function. 

a : 
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Q . II. DOMINANT SMALL k SINGULARITY 

>'. 

' ^ • In order to analyze the large w-behavior of the distribution function, it is convenient to use the Fourier transform 
' '/'(fc, of /(w, i). It is the generating function of the moments If /(w, t) has a tail ~ then the moments 

^ : with n > a are divergent, and so is the n-th derivative of the generating function at A: = 0, i.e. (j){k,t) is singular at 
■ - - fc = 0. Suppose the dominant small-fc singularity of 4>{k,t) is ~ where a is different from an even integer (even 
powers of k represent contributions that are regular at small k), then the inverse Fourier transform scales as l/|u|""'"'^ 
at large v. 

By applying Bobylev's Fourier transform method [Tl|jl^ , we obtain a nonlinear equation for 0(fc, t), and we deter- 
mine its dominant small-/c singularity. In doing so we have derived a transcendental equation for the exponent a in 
the power law tail ~ l/|c|°"'"'^ of the scaling solutions f{v/vo) of the Boltzmann equation for Maxwell molecules with 
in-elastic hard sphere interactions in arbitrary dimensions. 

The Boltzmann equation for the d-dimensional in-elastic Maxwell model reads. 
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dtfivut) = / dv2 [i/K*)/(«r) - /(«l)/(«2)] . (1) 



Here /^(- • •) = {1/fld) J dn{- ■ •) is an average over a d-dimensional solid angle where fid = 2n'^/^ /r{^d). The velocities 
V** with i,j = {1, 2|} denote the d-dimensional restituting velocities, and v* the corresponding posi-coZHsion velocities. 
They are defined as, 

V** = v,- i(l + -nn 
V* ^ Vi - ^{l + a)vij ■ nn, (2) 

with Vij — Vi — Vj, and n is a unit vector along the line of centers of the interacting particles. In one-dimension, 
the tensorial product nn can be replaced by 1. From the normalization of / it follows that the loss term reduces to 
—f(vi,t), i.e. the collision frequency is unity, and the dimension-less time t counts the average number of collisions 
per particle. 

We first illustrate the method for the one-dimensional case. Fourier transformation of the Boltzmann equation yield 
then, 

dt(b{k, t) = 4>{pk, t)4>{{\ - p)k, t) - <^(fc, t), (3) 

where we have used that 0(0, i) = 1 and p — \{1 + a). The equation for the scaling solution, 4>{k,t) — ^{vo{t)k) 
simplifies to, 

- -fkd^{k) /dk + $(fc) = $(pfc)$((l - p)k), (4) 

where the exponent 7 in VQ{t) ~ vq{0) exp[— 7t] is still to be determined. 

The requirement that the total energy be finite, imposes the the lower bound a > 2 on the exponent. We therefore 
make the ansatz that the dominant small-fc singularity has the form, 

$(fc) = l-i((fc.c)2)+^|A:|°. (5) 

Inserting this in (^, and equating the coefficient of equal powers of k yields the equation, 

1 -p" - (1 -p)" 



p(l -p) 



(6) 



The smallest root of this equation, satisfying a > 2, is a ~ 3, and A is left undetermined. Consequently the scaling 
solution has a power law tail, /(c) ~ l/c"*. 

The same method can be applied to the rf-dimensional case. Application of Bobylev's Fourier transform method to 
the Boltzmann equation for this case |l3,0] yields the transformed equation. 



dtmt)= / 0(fc+,t)0(fc_,t)-0(fc,t), (7) 

J n 

where the n-average is defined below (1), and 

fc+ — pk ■ nn = p^k'^(k ■ n)^ 

k_=k-k+ \k^\^ ^ k^[l ~ q(k ■ nf], (8) 

where q = p(2 ~ p) is a positive number. We proceed in the same way as in the one-dimensional case, and obtain the 
equation for the scaling solution. 



--^kd^{k)/dk + <^{k) ^ / $(fc+)$(fc_). (9) 

J n 

Inserting the ansatz (5) into (8), and equating the coefficients of equal powers of yields, 

a7(ifc • cr) = /(|fc • cr - ifc+ • cr - \k- ■ cd, (lo) 



for s = 2, a. In order to carry out the angular c-average in {\q ■ c|°) with q = {k, k+, fc_} we choose q as polar axis, 

and denote q-c = qcq-c — qc cosO, then {\q ■ c|") = \q\"- (Ic]"") Ki'^\ where K^'^'' is the average of | cos^l" over a 
d-dimensional solid angle, which equals 
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K(f^ = r(i(a + 1)) r(id)/r(i(a + d)) r(i), (ii) 

where q — p{2 — p). Finally we carry out the angular n-averages using (8), and obtain, 

|fc+|° ^ k^-p^Ki^^ 

= r /[I - q{k ■ nf]"/^ = k''L^^\q). (12) 



Insertion of these results in (10) for a=2, yields 

7=ip(l-rt = ij(l-«2). (13) 

For the exponent a, featuring in the power law tail of the scaling function /(c) ~ l/c°+'*, we obtain the transcendental 
equation, 

l-p'^K^''^-L^^\q) 
The two most interesting cases are d = 2, 3, where 

I.tt/2 



-flcos^ 611 



Jo 

L'i\q)= f dx[l-qx^T'\ (15) 

"'0 

and one can verify that a = 2 is also a solution of (14). We look for the smallest solution a{a) of this transcendental 
equation with a > 2. The numerical solutions for d = 2, 3 are shown in Figure]^ as a function of a. If p = i(l + a) t 1 
the root a{a) moves to oo, as it should, which is consistent with a Maxwellian tail distribution for the elastic case. 
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FIG. 1. Solution of Eq.(p^) as a function of a for 2 diemnsions (left panel) and 3 dimensions (right panel). The solution 
diverges as a ^ 1 (elastic limit), because / becomes Maxwellian in this limit. We note that a = 2 always satisfies Eq.(|l^, and 
the solution shown here is the one different from a = 2. 



The simulations in Ref. |y] of the two-dimensional Maxwell - Boltzmann equation show for the exponent 0(0; = 0) + 
2 ~ 5, where our analytical method predicts a{a = 0) + 2 ~ 6.2. In fact, closer inspection of their two-dimensional 
scaling plot shows that the slope of their log-log plot of /(c) versus c increases at larger velocities, approaching the 
exact prediction of the Boltzmann equation. However, at these large velocities the statistical errors in their simulations 
are too large to make a quantitative comparison for larger a-values. 
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